Reply to Comment of Gazdzicki and Heinz on Strangeness 

Enhancement in p + A and S + A * 

M. Gyulassy 1 , V. Topor Pop 1 - 2 , and X. N. Wang 3 
1 Department of Physics, Columbia University, New York, N.Y. 10027 
2 Institute for Space Sciences,P.O.Box MG - 6, Bucharest, Romania 
\ 3 Nuclear Science Division, Lawrence Berkeley Laboratory, 

ON 

Q\ ; University of California, Berkeley, CA 94720 

S3 

The Comment of Gazdzicki and Heinz |p| is flawed because their as- 
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sumed baryon stopping power in pA is inconsistent with data and because 
they ignored half the analysis based on the VENUS model. The Com- 
ment continues the misleading presentation of strangeness enhancement by 
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focusing on ratios of integrated yields. Those ratios discard essential exper- 
ts 

Q\ ■ imental information on the rapidity dependence of produced A and obscure 
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discrepancies between different data sets. Our conclusion remains that the 
NA35 minimum bias data on p+S — ► A+X indicate an anomalous enhance- 
ment of central rapidity strangeness in few nucleon reactions that points to 



^ . non-equilibrium dynamics as responsible for strangeness enhancement in 

H ; 

nuclear reactions. 

The Comment |l[ of Gazdzicki and Heinz addresses our recent analysis |§, where we 
concluded that if the NA35 data on p + S and S + S are correct, then strangeness 
enhancement must be due to new non-equilibrium multiparticle production dynamics. 
We further concluded that it is misleading to analyze strangeness enhancement in terms 
of ratios of integrated yields such as E s = (A + K + K) / (it) as advocated in Refs. 
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||.[4|]. Such ratios discard essential experimental information on the rapidity distribution 
of produced strange particles which contradict the conclusions drawn from those ratios 
alone. Such ratios also obscure striking discrepancies between different experiments which 
indicate that additional experiments, especially with pA, will be needed to understand 
the strangeness enhancement phenomenon. 

The focus of our work was the anomalous rapidity distribution of A reported by NA35 
|| for minimum bias p + S reactions. Our analysis used two different microscopic models, 
HIJING || and VENUS ||, to quantify the differences between p + p, p + A and A + A 
reactions. Gazdzicki and Heinz QXJ] missed entirely the VENUS half of our analysis and 
argued incorrectly and without calculation on the role of baryon stopping power to explain 
the pA data. 

In this reply, we contradict two key aspects of the argument presented in [|TJ: (1) their 
assertion that there is no strangeness enhancement in p + S reactions and (2) their asser- 
tion that our analysis should be dismissed because one of the models used, HIJING, 
underestimates the baryon stopping power of nuclei and overestimates the rapidity density 
of A's in the fragmentation regions of p + p. Other points raised in |IJ in connection with 
error bars on integrated yields and the pt acceptance cuts are irrelevant to our conclusions 
and will not be addressed here. 

First we review the NA35 data on minimum bias p + S || that contradict assertion 
(1). Then we review the known baryon stopping power from p + A — > p + X [^] and 
discuss the significance of the VENUS model simulations to contradict assertion (2). 

The primary motivation for our work was the very unusual rapidity distribution 
of A's reported by NA35 for the minimum bias p + S reactions. The central (y ~ 3) 
rapidity density (0.06 ±0.01) of A's produced in p + S was reported to be four times more 
than that (0.016 ± 0.0005) in p + p || (Fig la,b [Q). This is very surprising since the 
cross section for the selected p + S (n c h > 5) events was 470 mb, which corresponds to 
a minimum bias and not central trigger. Taking the diffuse nuclear surface into account 
that trigger only constrains events to have impact parameters less than 5 fm. For this 
weak trigger, the mean number of interacting S nucleons is only v « 2.2. The p + S 
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reaction therefore only tests the difference between multiparticle production dynamics in 
two nucleon [y = 1) and few nucleon {y — 2 — 3) reactions. 

/,From extensive p + A — > p + X systematics [[!]] the average baryon rapidity shift in 
p + A reactions grows slowly as Ay b ~ l + l)/3. lnp + S therefore the leading baryon 
rapidity shift is only a half a unit greater than in p + p. The qualitative discussion on 
baryon stopping in fl[ is wrong. Stopping power cannot account for the strong suppression 
of A production for y > 4.5 in p + S. Also the number of collisions v is too small 
to account for the factor of four enhancement of the central rapidity density. In the 
target fragmentation region (y « 1) there is also a factor of four enhancement of the A 
density in pS relative to pp. If the p + S data are correct, then strangeness enhancement 
already occurs in few nucleon processes and therefore must be due to new non-equilibrium 
dynamics. 

In stating the conclusion, we carefully pointed out however the fact that the NA35 
data on p + 5* differ substantially from earlier data on pAr by NA5 || and preliminary 
NA36 [|TIJ] pPb data. The earlier NA5 analysis [[| showed that both the Dual Parton 
and Lund models could account easily for the factor of two enhancement of the central 
A rapidity densities in p + At and similar enhancement in p + Xe and p + Xe. However 
the NA5 p + Ar — > A + X central density is a factor of two lower than found in p + S by 
NA35. The preliminary NA36 data also suggest that the mid-rapidity A density in p + Pb 
may be a factor of two lower than in p + S of NA35. Furthermore, the p + A — > A + X 



systematics of Ref. [T1J show that the enhancement of A production in p + A increases 
linearly with the number of secondary collisions. Therefore, strangeness enhancement 
builds up gradually in pA according to those other data, and is anomalously enhanced in 
the central region according to NA35 data. The ratios of integrated yields however hides 
these discrepancies completely, and that alone should be enough to disqualify such ratios 
as useful observables to study strangeness production in nuclear reactions. 

In we assumed the validity of the NA35 p + S data since otherwise the validity 
of the even more spectacular central S + S data would also have to be questioned (the 
NA35,NA36 discrepancy in central SS also remains unresolved). If, on the other hand, 



data on strangeness production eventually converge toward the lower NA5, NA36 range, 
then indeed our conclusions would have to change to a more pessimistic one - namely 
- conventional hadron cascade models could then account for most of the strangeness 
production (see e.g. Fig 13-15 jnj). In that case there would of course be no need to 
debate further this topic. 

Finally, we contradict assertion (2) of the comment Q. The overestimate of the 
fragmentation region A rapidity density with HIJING in p + p in spite of the accurate 
agreement with the mid-rapidity density is used to dismiss our model analysis. However, 
the fact that unlike in p + p both HIJING and VENUS fit the target fragmentation 
(y = 1) peak p + S in Fig lb |2| is itself direct evidence for strangeness enhancement. We 
note that the enhancement of target fragmentation region is consistent with NA5 and 



Ref. |TTJ and is not controversial. However, the strong suppression of A in the projectile 
fragmentation region (y > 4.5) in p + S is not found in either model. The authors of |IJ 
missed entirely the significance of the VENUS model calculation which reproduces well 
the p + A — > p + X stopping power measurements |J . In fact VENUS reproduces well 
the NA35 central SS — > pX distribution. HIJING is too strongly peaked about the mean 
rapidity loss while VENUS distributes baryons more broadly about that mean in accord 
with data. The small shift of the VENUS curves in the fragmentation regions (Fig. lb 
@) relative to HIJING show that correct baryon stopping cannot account for the strong 
suppression of A's with y > 4.5 in p + S. These calculations prove that the qualitative 
arguments used by the authors to account for the enhanced central rapidity A's three 
units of rapidity away from the fragmentation regions due to baryon stopping are wrong. 
Furthermore, as we emphasized in the agreement of the central rapidity density with 
pS data by VENUS is itself another independent proof that the central rapidity density 
in that reaction is anomalous since VENUS overestimates the p + p —>■ A + X central 
density by a factor of two. 

The analysis of ratios of integrated yields based on simplistic fireball models and con- 
cepts H,[l2| is misleading since the glaring discrepancies between the assumed differential 
distributions and data (and between data sets) can be well hidden from scrutiny. There 



can be no justification for throwing away valuable experimental information on differ- 
ential yields. Conclusions about the necessity of quark-gluon plasma production based 
on production ratio systematics alone are thus not well founded. We have shown in 
for example that at least one dynamical model (VENUS) can account for the anomalous 
central SS Lambda production distributions. Contrary to the assertion in JIJ strangeness 
enhancement does occur in pA at least in the target fragmentation region [jllj and pos- 
sibly in the central region as well |§. However, better data on pA will be required to 
resolve experimental discrepancies and to understand the non-equilibrium dynamical pro- 
cesses responsible for that enhancement. The search for unambiguous signatures of a new 
states of matter requires untangling complex and as yet poorly understood multiparticle 
dynamical effects that can forge naive plasma signatures. Only the detailed systematics 
and correlations between differential observables will be useful in that search. 
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